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ABSTRACT 

A new model is presented which explains the origin of the broad emission 
lines observed in the LINER/Seyfert nucleus of M81 in terms of a steady state 
spherically symmetric inflow, amounting to ~1 x 10~^ Mq/jt, which is sufficient 
to explain the luminosity of the AGN. The emitting volume has an outer radius of 
~1 pc, making it the largest broad line region yet to be measured, and it contains 
a total mass of ~ 5 x 10~^ Mq of dense, ~ 10^ cm~^, ionized gas, leading to a 
very low filling factor of ~ 5 x 10^^. The fact that the BLR in M81 is so large 
may explain why the AGN is unable to sustain the ionization seen there. Thus, 
the AGN in M81 is not simply a scaled down quasar. 

Subject headings: galaxies: Seyfert, galaxies: individual(M81, NGC 3031), quasars: 
emission lines 



Introduction 



Broad emission lines (ISeyfert Ill943l ) are a defining characteristic of active galactic nuclei 
(AGN). They occur in AGNs spanning seven decades in luminosity, from the lowest lumi- 
nosity Seyferts to the highest luminosity quasars. A variety of models have been proposed 
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to explain the origin of the broad emission hnes, but they all have one feature in common; 
that the lines are formed in the dense, ~ 10^ to 10^ cm~^, and fast moving, ~ 10'* km/s, 
gas, close to a massive black hole (MBH) (iReed 119771 ) . The models differ only in how the 



gas achieves it's high velocity, be it through inflow, outflow, spinning in a disk or attached 
to stars orbiting close to the AGN. Identifying the mechanisms that produce broad emission 
lines is of paramount importance as only then will we finally be able to understand, and 
exploit, the physics of the small and as yet, spatially unresolved broad emission line regions 
(BELRs) in galaxies. 

The main purpose of this paper is to draw attention to an important subset of AGNs 
that exhibit fiat-topped broad emission line profiles which may be utilized to determine a 
size for the broad line region (BL R) that her e tofore has previously o nly been possible using 
reverberation mapping techniques (iPeterson l Jl993h . iPeterson I (120011 ) ). We demonstrate the 
utility of this new technique by interpreting the fiat-topped broad Ha emission hne profile 
of the LINER/Seyfert nucleus of M81. 



At a distance of 3.6 Mpc (IFreedman et al. 1120011) M81 contai n s the nearest and best 



resolv ed low luminosity AGN (jPeimbert fc Torres-Peimbert I (1198 ll ). iFilippenko fc Sargent 



(119881 )). The almost face-on aspect of the host galaxy allows a clear and unobscured view of 
the nu c leus which has been studied across t he electromag r ietic s p ectrum (iPeimbert fc Torres-Peimbert 
Ji98i[ ). IHo. Filippenko fc Sargent Nl996h . lBower et al. I (fl996[ ). lBietenholz. Bartel fc Rupin 



(120001 ) ■ iLa Parola et al. I (120041 )1 M81 harbors a black hole (BH) with a mass, M^^ 



X 10^ M(7^, deduced fr om a two dimensional velocity map obtained with STIS on the HST 



( jPevereux et al.ll2003l ). The AGN is identified wit h a UV point source at the nucleu s which 
has a luminosity at 1500 A, L1500 = 5 x 10^ Lf. 
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Our paper is based on the premise that the broad emission lines (BELs) are produced in 
a kinematically distinct and hence physically distinct region centered on the BH. Therefore, 
the broad emission line region (BELR) has a well defined inner radius and outer radius. 
Since the gravitational field strength is already known in M81, the relationship between 
velocity and radius may be established, given a kinematic model for the BELR gas. Thus, 
in principle, one can determine the inner and outer radii of the BELR in M81 by modeling 
the broad emission line profiles. 

The layout of our paper is as follows. In section 3, we evaluate the origin of the BELs 
seen in M81 in the context of viable models for the production of BELRs in AGNs, namely 
inflow, outflow, rotating gas disks and the atmospheres of stars illuminated by the AGN. In 
section 4, we discuss the viable models in the context of various constraints and conclude, 
in section 5, that inflow provides the best explanation for the BELR in M81. We begin. 
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however, with a review of the broad Ha emission hne observed in M81. 



The Broad Ha Emission Line Profile in M81 



The br oad Ha emission hne , measured by STIS aboard the HST, has been presented 
previously (jPevereux et al.ll2003l ) and is shown in Figure 1 for convenience. As the figure 
illustrates, the profile exhibits a rather unusual 'fiat-top', that is rarely seen in AGNs. Such 
profiles may be produced by radiation from a spherica. 



l y sym metric shell of gas in radial 
motion as noted previously by iCapriotti. Foltz fc Bvard I (|l980r). Single pe aked emission line 
profiles may also be produced by accretion disks (jchen &: Halpern 1989 ) and bipolar flows 



fIZheng. Binette k Sulentic Ill990l ). 



As reported in iDevereux et al.l ( 20031 ). the entire broad Ha line is redshifted by ~ 7 A 
(~ 300 km/s) with respect to the narrow Ha line. The broad line is very symmetrical and so 
the redshift can not be attributed to skewness ii i the line shape. Similar red shifts have been 
noted for broad emission lines in other AGNs (lEracleous fc Halpern II2003I ). We presently 
have no explanation for this phenomenon, and we overlook it for now, with the consequence 
that the velocities reported in this paper are relative to the center of the respective emission 
lines. 

Slight irregularities on the redward side of the proflle, near the peak, coincide with 
the narrow components of the 6563 A Ha line and the brighter of the two [Nil] lines, 
each of which had to be subtr a cted to isolate the broad Ha component, as described in 
more detail by IDevereux et al.l (120031 ). The blue side of the proflle, however, is less sus- 
ceptible to the subtraction procedure because it coincides with the fainter [Nil] line. In- 
terestingly, the H/? lin e is sim ilar in shape to the Ha emission line as noted previously by 
Filippenko &: Sargent I (119881 ) and lBower et al. I (119961 ). The distinctive 'flat-top' is thus well 
establish ed for M81. Other broad emissio n lines are observed in M81, as described previ- 
ously by IHo. Filippenko &: Sargent I (Il996l ). but the Balmer lines are the brightest, they are 
presently the best resolved spectrally and hence the easiest to model. 

The broad Balmer lines in M81 are apparently variable. Faint, widely sep arated 'double 
peaks ' did appear on the wings of the Balmer lines as reported previously by iBower et al. 
(119961 ) ■ but they did not appear in the more recent STIS spectrum. The transient features 
are consistent with the disruption of a star at the tidal radius of the MBH in M81, but the 
broad line itself remains to be explained. 
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3. Broad Line Region Models 



A large body of work exists addressing the o rigin of broad lines in lumi nous AGNs, but 
the problem of the fo rmation and confinem ent (IMathews &: Ferland Ill987l). the excitation 
( iKwan &: KrolikI Il98ll ) and the kinematics (iRobinson. Perez. &: Binettd llQQOl ) of the broad 
l ine gas remains unsolved. Interestingly, broad lines also occur in low luminosity AGNs 
( iFilippenko fc Sargent Ill985l ) which presents the opportunity to examine the phenomenon 
from a different perspective. Furthermore, rather than attempting to explain the observed 
properties of all low luminosity AGNs, we approach the problem from the point of view of 
examining one well studied example, namely M81. 



3.1. Outflow Model 



Ev idence for ou t flow is seen in the spectra of broad absorption line quasars (e.g. [Hazard et al 



(119841 ). lArav et al. I (119991 )). However, the outflow model can be ruled out for M81 because 
the diminutive luminosity generated by the AGN is simply unable to provide sufficient ra- 
diation pressure to overcome the gravitational force of the MBH as demonstrated in the 
following. 



The condition for a radiatively driven wind is given by. 



KL/AncG > M 



where M and L are the mass and luminosity respectively, interior to a radius r, c is the 
speed of light, and G the gravitational constant. For M81, this condition is not satisfied by 
six orders of magnitude, if one adopts the Thompson scattering opacity k = 0.4 cm^/g for 
a pure hydrogen gas, M^h = 7 x 10'' Mq, and a combined UV and X-ray luminosity of 5 
X 10^ L0. The disparity is simply too large to overcome even by invoking line opacities, as 
line d riven winds are y i able only for objects that radiate close to t he Eddington luminosity 



one sided radio jet ( 



limit (IKing fc Pounds I (120031). IShlosman. Vitello fc Shaviv I (Il985h l A short, 1.7 x 10"^ 



Bietenholz. Bartel fc Rupin Il2000l ) is seen to emerge from the nucleus 



of M81, but this is an unlikely source for the observed broad line emission as that outfiow 
consists of a relativistic plasma. 
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3.2. Stellar Atmospheres Illuminated by the Central AGN Model 



We are also able to rule out another contending theory for the production of broad 
emission lines in M81, namely ion ization of the extended mas s loss envelopes of red gi a nt sta rs 
orbiting close to the AGN (eg. IScoville fc Norman I (119881 ) . [Alexander fc Netzer I (jl997l )). 
Stellar winds provide a natural explanation for the confinement and replenishment of the s o 
called 'broad line clouds', which is otherwise a major problem (IMathews fc Ferland 1 119871 ) . 
However, such a model is unlikely to be the explanation for the BELR in M81, because the 
intensity of ionizing photons at the tidal radius is insufficient to penetrate the denser, ~ 10^ 
cm~^, region of the mass loss wind, as shown in the following. 

The tidal radius, Kudai, for stars in synchronous rotation around a BH is given by 



Rtidal — Rstar[M BH / M 



star\ 



11/3 



(2) 



For a 1 M© star of radius 10^^ cm, in the vicinity of a 7 x 10^ M© BH, the tidal radius 



corresponds to 4 x 10^^ cm 



Following IScoville &: Norman I (119881 ) , and assuming plane parallel geometry, the ion- 
izing photon flux at the tidal radius is related to the penetration depth, d, in the inverted 
Stromgren sphere of the stellar wind by 



[L/{A7TRudai^hiy)]d 



V 



(3) 



where olb = 2.6 x 10 cm ^ s ^, is the Case B recombination coefficient, and hv is 
the energy of an ionizing photon. 

For a power law of spectral index a. 



(4) 



the right hand side of equation 3 integrates to 



/■oo 

/ aBuj^dr = LoU'^[v~°' - u~'^^]/{ha4nRtidai^) 
Jd 



(5) 



The hydrogen number density, n^, is related to the stellar mass loss rate, M, and the 
wind velocity, vw, via the equation of continuity; 
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(6) 



Combining equations 3 and 6, and evaluating the left hand integral, one finds, 



1/3 



(7) 



Subst ituting equation 7 into equation 6, and adopting a = 2, hug = 13.6 eV and hz/^ax 
= 600 eV (IHo. Filippenko fc Sargent Ill996l ). one obtains for the number density in the wind, 
n, at the penetration depth, d, of the ionizing photons; 



n{d) = Q3xlO^[Li5oo/W^LQ]'^^^[Rudai/W^^cm]-^^^[vw/lOkm/s]^/^[M/10-^MQ/yr]^^^^ 



Substituting representative values for the wind velocity Vw = 10 km/s, and stellar 
mass loss rates of M = 10"^ Mq/jt JScoville fc Norman I Hossh and M = 10'^ Mq/jt 



( lAlexander &: Netzer 1 119971 ) . one finds for M81 that representative densities in the wind 
at the penetration depth of the ionizing photons are 30 x 10^ cm~^ and 63 x 10^ cm~^, 
respectively. Thus, the density where most of the Ha emission is produced is expected to 
be significantly lower than the ~ 10^ cm~^ required for the bro ad line region of M81 to not 



produce broad forbidden lines (IHo. Filippenko fc Sargent Ill996l ). Furthermore, the densities 
calculated above are for stars as close to the UV source as the tidal radius permits, and 
assume no absorption of ionizing photons inside that volume. Obviously, the penetration 
depths and densities will be lower for the majority of stars which are further away. Plus, 
it is quite likely that there is gas inside the tidal radius that would attenuate the ionizing 
fiux, thereby exacerbating the problem. An additional argument against stars is that, for 
M81, the emitting envelopes would occupy such a narrow range of orbital velocities that they 
would be unable to reproduce the observed broad line profile. We therefore conclude that 
the BELR in M81 is not produced by the extended mass loss envelopes of red giant stars 
orbiting close to the AGN. 



3.3. Accretion Disk Model 



The next contender for the origin of the b road emission lines is a rotating 'accretion 
disk'. While there has been a novel proposal (IMurry fc Chiang 1 119971 ) to generate a sin- 
gle peaked profile using a disk wind, we are reluctant to adopt it as there is no a priori 
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evidence that there is an outflowing wind in M81. Instead, the m odel used here involves 
a relativistic accretion disk developed by I Chen &: Halpern I (Il989l ). and adopted later by 
Eracleous fc Halpern I (120011 ) to explain the single peaked broad line profile in the LINER 
nucleus NGC 3065. An illustrative fit to the broad Ha emission line in M81, using the 
relativistic disk model, is presented in Figure 2. 

The model invokes five free parameters, a dimensionless inner radius, (i, and outer 
radius, an inclination angle, i, an emissivity law, q, and a velocity dispersion for the 
gas, a. As described previously by lEracleous fc Halpern I (120011 ) . the two emission peaks; 
a distinguishing characteristic of disk models, may be merged into a single peak, and the 
profile wings made symmetric, if the disk has a large outer radius and is inclined to the line of 
sight. The model also incorporates a 300 km/s velocity dispersion to widen the two velocity 
peaks so that they merge into a reasonable representation of the flat-top profile seen in M81. 
Additionally, the model incorporates a q=2 emissivity law which accentuates emission from 
the inner regions of the disk. For a BH with a mass the same as measured for M81, the outer 
radius of the disk model, illustrated in Figure 2, is large, corresponding to ~ 0.8 pc. 

The success of the model in fitting the profile is dependen t on the disk iriclinat ion. If 
the disk is nearly face-on, as observed for the narrow line gas ( iDevereux et al.ll2003l ). then 
the inner radius of the disk has to be small to achieve the necessary range of velocities 
at zero intensity. Consequently, gravitational and transverse redshifts become discernable, 
with the result that the red wing becomes asymmetric, somewhat extended and unlike the 
observed profile. However, a disk inclination of 50 degrees is able to reproduce the shape of 
the M81 profile quite well as Figure 2 illustrat es. Such a high inclination for the inner disk 
is supported by observations of the radio jet (IBietenholz. Bartel fc Rupin 1120001 ). Thus, if 
the broad emission line is produced by an accretion disk, then the unresolved inner part of 
the disk is considerably more inclined than the larger scale disk seen in ground based and 
space based images. 



3.4. Evidence for Inflowing Gas In M81 

As alluded to in the introduction, a steady state spherical inflow is also able to reproduce 
the broad Ha emission line profile observed for M81, and with the minimum of assumptions. 
The relationship between velocity and radius is determined by the mass distribution, M(r), 
to be 



v{r) = 



(9) 
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where M(r) includes a point mass, representing t he black hole, erabedd ed in the center 
of an extended star cluster, as described previously by lDevereux et al.l (120031 ) . The observed, 
radial component, of the velocity for each particle, i, is given by 



Vi = v{r)cosQ (10) 

where G is the angle between the radial velocity vector and the line of sight. 

A broad line profile, $(w), with the characteristics of the one observed in M81 may 
be produced by generating a histogram of velocities for a system of particles distributed 
randomly within a spherical volume bounded by an inner radius. Tinner, an outer radius, 
i^outer and a cloud number density distribution, N(r), such that, 

"^outer 

For a steady state flow, where M ^ f(r), the cloud number density distribution is 
determined by mass conservation to be. 



N{r) oc r~^/2 



(12) 



As described previously by ICapriotti. Foltz fc Byard I (Il980l ). the velocity width of the 
'flat-top' is determined by the outer radius, and the width at zero intensity is determined 
by the inner radius. Thus, given a density distribution and a velocity law, one can use the 
profile shape to determine the physical size of the emitting region. Values that reproduce 
the best representation of the broad Ha emission line profile, shown in Figure 1, correspond 
to an outer radius, router = 1 pc, and an inner radius, r^ner = 0.2 pc. 

There are other broad lines seen in M81, but they are fainter and were acquired with 
lower spectral resolution. However, two of the better examples appear in the UV; Lya at 1216 
A and Mg H at 2800 A. They ar e reproduced in Figures 3 and 4 respectively. According to 
Ho. Filippenko fc Sargent I (Il996l ) , the central depression seen in the Lya line is likely caused 
by interstellar absorption and the bright central narrow line is geocoronal emission. They 
suggest that interstellar absorption is also responsible for the structure seen in the core of 
the Mg H line. Nevertheless, the wings of the Lya and Mg II profiles can be reproduced 
using the same inflow model and yield slightly different sizes for the emitting region. For 
example, the Lya emission hne requires router = 0.7 pc, and an inner radius, r^ner = 0.1 pc. 
The Mg II emission line requires router = 0.9 pc, and an inner radius, runner = 0.35 pc. These 
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illustrative fits are shown in Figures 3 and 4. Overall, the sizes determined for the BELRs 
are similar. We note that the line width at the intensity of the continuum sets an upper 
limit on the inner radius. This is because the line may well persist below the continuum 
level, even though it can not be directly measured. 



Discussion 



Evidently, two plausible models; a rotating accretion disk and a spherically symmetric 
inflow, are able to reproduce the shape of the broad Ha emission line profile observed for 
M81. In either case, the velocity law dictates that the emitting region is large; ~ 0.7 to 1 pc 
in radius. In the context of the accretion disk model, the broad Ha emission line presumably 
represen ts the continuation of the larger s cale Ha disk, seen in ground based and space based 
images (jPevereux. Ford fc Jacoby 1 119971 ) , as it extends down into the unresolved nuclear 
region. On the other hand, the interpretation of the broad Ha emission line profile in terms 
of a spherically symmetric inflow would represent a new phenomenon for M81. 



4.1. Broad Line Region Size 



Perhaps the most surprising result to have emerged from our analysis is the large outer 
radius, ~ 1 pc, inferred for BELR of M81, regardless of whether the broad emission line 
is produced by an accretion disk or an inflow. With an angular diameter of 0.1", the BLR 
of M81 is on the brink of being resolved with HST. The large outer diameter inferred for 
the BLR may explain why a 'flat-top' proflle is observed in M81, and rarely in other AGNs. 
Unless the BELR is resolved, the line profiles will be contaminated by the slower moving gas 
surrounding the BELR, causing the line profiles to be peaked rather than fiat. 

The size we infer for the BELR causes M81 to not conform to the correlation between 



BLR size and U V luminosity established for higher lu r ninosi ty AGNs ( iPeterson I (120011 ). 
Peterson I (119931 )). However, as noted by iKaspi. et al. I (120051 ) . the correlation appears to 
break down for low luminosity AGNs, with Ljjy < 10^ Lq, which would include M81. But 
the large size that we infer for the BELR of M81 does make it the largest ever measured. 
On the other hand, there are good reasons to believe the sizes of broad line regions may 
have been underestimated, previously, using the reverberation m apping technique, w hich is 



most sensitive to na e asurin g the inner radius, closest to the AGN (lEdelson fc Krolik I (119881 ). 
Robinson fc Perez I Jl99o[ )l 
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4.2. Broad Line Region Ionization 



The size we i nfer for the BELR in M81 is ord ers of magnitude larger than the range 
of values cited by IHo. Filippenko fc Sargent I (Il996l ) but their estimate is based on the ex- 
citation parameter which links the excit ation of the broad emission lines to th e central UV 
source . However, as noted previously by iHo. Filippenko fc Sargent I Jl996h andlMaoz et al 



( 1l998l ) the nuclear UV source in M81 produces insufficient ionizing photons to explain the 
luminosity of the broad Ha line, by a factor of 11 according to our calculations. We find 
that the the central AGN produces 9 x lO''^ ionizing ph/s, integrated between 13.6 and 600 
eV, assuming a po wer law, whereas excitation of th e broad Ha emission line requires 
1 X 10^^ ionizing ph/s. iHo. Filippenko fc Sargent I (119961 ) also note that the BLR does not 
respond to time variable X-ray emission. In our view, this provides further evidence that the 
central AGN is not responsible for ionizing the BLR gas. That the BELR in M81 is so large 
may explain why the central AGN is unable to sustain the ionization seen there. Thus, the 
low luminosity AGN in M81 is not simply a scaled down quasar. 

The origin of the excitation for the BELR is as enigmatic as the or igin of the excitation 

for the more extended narrow line region in M81, discussed previously by lDevereux. Ford fc Jacoby 
( 1l997l ). Most likely, the gas is excited, in situ, by shocks, or by hot, and as yet undetected, 
PAGB stars. In either case, the excitation parameter would then be completely unrelated 
to the central source, and hence would not yield a useful measure of the BELR size. 



4.3. Virial Black Hole Masses 



Previous attempts to estimate t h e BH mass in M81 using the s o called ^Virial Method^ 
( Peirnbert fc Torres- Peimbert I (ll98ll ). lFilippenko fc Sargent I (I1988[ ). IHo. Filippenko fc Sargent 



(119961) ) have consistently underestimated the kinematically determined mass (IDevereux et al. 
2003h by factors ranging from 20 to 200. For example, the recent formalism of I Greene fc Ho 



( 120051 ). which uses the FWHM and luminosity of the broad Ha emission line, underestimates 
the mass of the BH in M81 by a factor of 140. On the oth er hand, the BH in M81 d oes 
conform to the BH niass-bulge veloc i ty disp ersion correlation (iFerrarese fc Merritt Il2000l ) as 
noted previously by iDevereux et al.l (120031 ) . This dichotomy is regarded as further evidence 
that the BLR in M81 is very different from those studied in more luminous AGNs. 



4.4. The Mass of Ionized Gas in the BELR of M81 



The mass of emitting gas may be deduced from standard (Case B) recombination theory; 
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M, 



emitting 



(13) 



Using an effective recombination coefficient 

-3 



8.6 X 10^ 



-14 



cm s , assuming a 

constant average density n = 10^ cm"*^, and a lumin osity L (Hq,) = 3.7 x 10^ Lq based on 
the broad line ffux reported in Devereux et al. (j2003 ). one finds "^emitting = 0.05 Mq. This 
mass is implausibly small for an accretion disk, but can be more easily reconciled with an 
inflow as shown in the following. 



4.5. The Filling Factor and the Inflow Rate for the BELR of M81 

It is straight forward to calculate the filling factor, e, once the dimensions of the emitting 
region have been established. For a uniform density medium occupying a spherical volume 
of radius r, one finds 



?,L{H^)lA'Knla%ihuHy 



(14) 



Again, using an effective recombination coefficient oi^hL ~ ^-^ ^ 10 

-3 



14 



cm s , assuming 

a constant average gas density n = 10^ cm~^, and a lumi nosity L (Hq.) = 3.7 x 10^ based 
on the broad line flux reported in lDevereux et al.l (120031 ) . one flnds e ~ 5 x 10"^ for M81. 



Having established the dimensions of the emitting region and the flUing factor one can 
now calculate the mass inflow rate for the ionized gas, using the equation of continuity; 



M = eA'Kr'^v{r)nH{r)mH (15) 

The velocity at the inner radius of 0.2 pc is determined by the mass distribution to be 
1764 km/s. Setting the gas density in the flow to be 10^ cm~^ one obtains a mass inflow 
rate, M = 1 x 10~^ VIq/yi. If the gas density is higher, the mass inflow rate will decrease. 
Conversely, if only a fraction of the inflowing gas is ionized, then the total mass inflow rate 
will, of course, be higher. 



4.6. The Luminosity of the AGN 



A variety of models have been considered for the production of the broad Ha emission 
line in M81, but by far the most straight forward model, that involves the minimum of 
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assumptions, is the inflow model. Accretion onto a compact massive object has long been 



suspe cted as the origin for the luminosity of the AGN in M81 (jPeimbert &: Torres- Peimbert 



198ll ) and recognizing the shape of the broad line proflle may be the signature of an inflow 



permits us to corroborate this idea. According to the virial theorem, the maximum power, 
P, produced from the conversion of gravitational potential energy as material falls from a 
distance of Ipc, which for the purposes of this calculation is effectively at inflnity, onto the 
event horizon of a MBH, can be expressed as 



P = Mc^/A 



(16) 



A lower limit to the mass inflow rate that we infer from the ionized gas, 10~^ M0/yr, sets 
a lower limit of 4 x 10^ for the power, P, which, without explaining the details, is sufficient 
to account for both the UV and X-ray luminosity of the AGN with a conversion efficiency of 
power into radiative luminosity of < 14%. Interestin gly, the mass inflow rate tha t we in fer 
is also sufficient to sustain the radio jet according to lBietenholz. Bartel fc Rupin I (120001 ). 



4.7. The Origin of the Broad Line Gas 



If the broad emission lines are produced by an inflow, then one wonders how the incred- 
ibly small mass of gas in the BELR, measured to be 5 x 10~^ Mq, is distributed within the, 
by comparison, rather large emitting volume ~1 pc in radius. Oi ie can get an approximate 
estimate based on the 'noise' in the proflle, as flrst suggested by ICapriotti. Foltz &: Byard 
(fl981i) . The fact that the observed broad line proflle is smoother than the model proflle, 
which was generated using 75,000 points (see Figure 1) suggests that, on average, the mass 
of each 'broad line cloud' is < 10~^ M0. We surmise that repres e ntativ e masses for broad 
line clouds ~ 10^^ M©, cited previously by [Alexander &: Netzer I (119941 ). are not at all un- 
reasonable given that the observed proflle in M81 is, in fact, so smooth. 

It is not surprising, given the small mass of gas and the ver y low volume fllling factor, 
that the internal extinction is negligible to both the UV sou rce (IHo. Filippenko &: Sargent 
19961) and th e gas responsible for the broad Balmer lines (iFilippenko fc Sargent I Jl988h , 



Bower et al. 



( 19961 )). The fact that the Ha proflle is symmetric provides further evidence 
that the e xtinction to the BELR is i ncons equential; apparently a common feature among 
AGNs (eg. IColhn-Souffrin fc Lasota I Jl988h l 



We speculate that mass loss from a spherically symmetric distribution of late-type bulge 
stars may provide a more than adequate supply for the in-falling gas. Any gas for which the 
vector cross product 



rxp = 



(17) 



where r is the radius vector and p is the hnear momentum of the gas, will fall into the 
nucleus. Calculating how much gas satisfies this identity and furthermore, how it becomes 
so compressed, is beyond the scope of the present paper, but we note that only a small mass 
inflow rate, ~ 10~^ Mo/yr, of dense, 10^ cm~^, gas is required to explain the luminosity of 
the AGN. 



5. Conclusions 

We have demonstrated a new technique that has allowed us to determine the size of 
the broad hne region (BLR) in M81 by modehng the shape of the broad hne profile using a 
kinematic model for the broad emission line gas. Our principle conclusion is that the BLR 
is large, ~ 2 pc in diameter, regardless of whether the velocity field is represented by an 
accretion disk or a spherically symmetric inflow. 

The fact that the BLR in M81 is so large may explain why the AGN is unable to 
sustain the ionization seen there. We therefore conclude that the gas responsible for the 
broad emission lines seen in M81 must be ionized in-situ, either by shocks, hot stars, or a 
combination of the two. Thus, the AGN in M81 is not simply a scaled down quasar. 

How the gas becomes so compressed as to prevent the formation of broad forbidden lines 
remains a mystery. Nevertheless, if the gas density really is high, ~ 10^ cm~^, then the broad 
Ha emission hne observed in the LINER/Seyfert nucleus of M81 is most easily understood in 
terms of a steady state spherically symmetric inflow, amounting to > 1 x 10~^ Mq/jt, which 
is sufficient to sustain the luminosity of the AGN and the radio jet. We further speculate, 
based on the smoothness of the broad emission line profile, that the individual regions of 
dense ionized gas have masses < 10~^ Mq. Alternatively, if the gas density is significantly 
lower, then the total mass of gas increases, making the interpretation of the broad emission 
lines in terms of an accretion disk more amenable, but then explaining the absence of broad 
forbidden lines becomes a problem. Thus, the gas density remains pivotal in understanding 
the origin of the broad emission lines seen in M81. 
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Fig. 1. — Broad Ha emission line in M81. The observed profile (jPevereux et al 



2003' 



3) is shown in orange and the inflow model is shown in green. The inflow was 
generated using 75,000 randomly distributed points in free-fall from an outer 
radius of 1 pc (see text for details). The intensity of both proflles has been 
normalized to the intensity at zero velocity. 




Fig. 2. — Observed Broad Ha emission line in M81, shown in orange, with the 
disk model shown in green. The disk is 0.8 pc in radius and is inclined by 50 
degrees to the line of sight (see text for details). The intensity of both profiles 
has been normalized to the intensity at zero velocity. 
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Broad Lyg emissi on line in M81. 



The observed profile 
( IHo. Filippenko &: Sargent Ill996l ) is shown in orange. The model profile, shown 
in green, represents an infiow (see text for details). The intensity of the model 
profile has been normalized to the intensity at zero velocity. 
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Broad Mg II emis sion line in M81. The observed profile 



(IHo. Filippenko &: Sargent Ill996l ) is shown in orange. The model profile, shown 
in green, represents an infiow (see text for details). The intensity of both profiles 
has been normalized to the intensity at zero velocity. 



